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Abstract: In an extension of the minimal supersymmetric standard model (popularly 
known as the /Uz^SSM), three right handed neutrino superfields are introduced to solve the 
/i-problem and to accommodate the non-vanishing neutrino masses and mixing. Neutrino 
masses at the tree level are generated through i?— parity violation and seesaw mechanism. 
We have analyzed the full effect of one-loop contributions to the neutrino mass matrix. 
We show that the current three flavour global neutrino data can be accommodated in 
the ^i^SSM, for both the tree level and one-loop corrected analyses. We find that it is 
relatively easier to accommodate the normal hierarchical mass pattern compared to the 
inverted hierarchical or quasi-degenerate case, when one-loop corrections are included. 
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1. Introduction 



Despite its stupendous success in explaining elementary particle interactions, the celebrated 
standard model (SM) of particle physics suffers from some shortcomings, both theoretical 
and experimental. On the experimental side, explaining the masses and the mixing pattern 
of neutrinos is a task in which the SM is an apparent failure. The neutrino sector, therefore, 
is a natural testing ground for most proposals for going beyond the SM. 

Supersymmetry (SUSY) is a rather popular choice for new physics. The minimal super- 
symmetric version of the SM (MSSM) provides a natural solution to the so-called "gauge- 
hierarchy problem" through the introduction of superpartners of SM-particles. However, 
MSSM itself is not free of drawbacks. One of these is the so-called //-problem [1], which 
essentially means our lack of understanding as to why the higgsino mass parameter, a 
SUSY invariant quantity, has to be around the SUSY breaking scale. This problem can be 
solved in the next-to minimal supersymmetric version of the standard model (NMSSM), 
where, unlike in the MSSM, the fi-tevm becomes a derived quantity with the right order of 
magnitude. 

Unfortunately, neither the MSSM nor NMSSM by itself can explain the observed pat- 
tern of neutrino masses and mixing, on which definite guidelines have been set down by 
existing data [2,3]. The situation becomes different if one allows violation of the discrete 
symmetry known as i?-parity (Rp) [4-12], defined as Rp = (_i)^+3S+2S'^ where L{B) is 
the lepton(baryon) number and S is the spin, of a particle. Neutrino masses and mixing 
have been analyzed in these models under various assumptions, both at the tree level and 
by taking loop-induced effects into account. Neutrino mass generation in variants of Rp 
violating MSSM have been addressed in refs. [13-62]. Such extensive study has estab- 
lished i?p-violation to be as potent in neutrino mass generation as the well-known seesaw 
mechanism [63-67] which requires introduction of gauge singlet neutrino superfields. 

R-parity violation, in both the contexts of accelerator phenomenology and, for example, 
neutrino mass generation, have been extensively studied in various scenarios, especially in 
its L-violating incarnation. Thus one has so-called trilinear R-parity violation driven by the 
Xijk or A^^-^-type terms in the supcrpotcntial. In addition to new signals induced by three- 
body decays of the lightest neutralino, neutrino masses are generated through loop effects in 
such scenarios [13 20]. Then one can have R-parity broken by bilinear terms [21-58] of the 
type ejLjiJ2- One here notices remarkable features like (a) non-zero vacuum expectation 
values for sneutrinos, and (b) the mixing between neutrinos and neutralinos as well as 
charged leptons and charginos. Such a scenario generates one neutrino mass at the tree level 
while the other mass(es) need to be generated via loop effects [39-51]. The characteristic 
signal consists in final state with comparable numbers of muons and taus at high energy 
colhders [52-58]. 

Bilinear R-parity violation, is also linked with spontaneous L-breaking [10, 59-62] via 
a singlet sncutrino vacuum expectation value (VEV). Such an effect triggers terms of the 
type eiLiH2 in the superpotential, and also leaves as its footprint a Majoron which has its 
own experimental signature [68-70], being an additional source of missing energy. 

In the backdrop of such a rich phenomenology of relatively minimalistic R-parity break- 
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ing models, further embellishments on the minimal scenarios have also been studied, often 
with some specific goals. A proposal for neutrino mass generation together with a solution 
to the //-problem, with the same set of gauge-singlet right chiral neutrino superfields, has 
been advocated in [71]. This is popularly known as fiuSSM. Following this proposal, the 
scalar sector and the parameter space of /iz/SSM was studied in great detail in [72]. In 
/xz/SSM, the three generations of SM-neutrinos can acquire masses through a TeV-scale 
seesaw mechanism, with both neutralinos and heavy neutrinos participating in the pro- 
cess [73] . Issues of neutrino mass generation and the //-problem can also be found in some 
recent works [74-78]. 

A comprehensive analytical study of mass generation of light neutrinos in /tz/SSM, 
accompanied by necessary numerical analysis, has been discussed in ref. [73]. In this work, 
neutrino masses and mixing, consistent with the three flavour global neutrino data, are 
reproduced even with the simplistic choice of flavour diagonal neutrino Yukawa couplings 
(yj*). Decay modes of the lightest neutralino into two-body final states [Z^vi^ W^£) 
have also been considered for various compositions of the lightest neutralino. In addition, 
correlations between neutrino mixing angles, and ratio of the decay branching ratios into 
W=''-charged lepton are studied there as a possible test of this model at the Large Hadron 
Collider (LHC). 

Among other related studies of importance, neutrino mass generation and collider 
aspects of this model, with one and two generation(s) of the right handed neutrinos, have 
been addressed in ref. [79]. For one right handed neutrino, neutrino mass generation has 
been studied there upto one-loop level. Decays of the lightest neutralino into all possible 
final states are also studied in this reference. Constraints on complex vacuum expectation 
values (VEVs) from clcctroweak symmetry breaking (EWSB), and its consequence on the 
neutrino sector are studied in ref. [80], where issues concerning spontaneous CP-violation 
are also addressed. The role of gravitino as a dark matter candidate in this model was 
studied in ref. [81]. This paper also highlighted the prospects of detecting gamma rays 
from decaying gravitinos. For an overview of various other aspects of /ii/SSM, see the 
recent review [82]. 

In this work, we study in detail the effect of radiative corrections, upto one-loop, to 
the neutrino masses and mixing, consistent with the three flavour global data [2,3]. As 
mentioned before, a similar study, but with just one generation of right handed neutrino, 
was carried out in [79]. However, a full study addressing both neutrino masses and the 
bilarge mixing pattern, with a complete set of one-loop corrections with all three genera- 
tions of left and right handed neutrinos, has so far been lacking. This is exactly what we 
attempt in the present work. We also perform a systematic study to identify the crucial 
parameters of the model, which control the tree level or the one-loop dominance in the 
neutrino sector. 

As shown in Ref. [73,79], a very attractive feature of this model is that the ratios of 
certain decay branching ratios show very nice correlation with the neutrino mixing angles. 
This is very similar to bilinear R-parity violating models [52,53,56]. Nevertheless, one 
should note certain differences in these two cases. In /xi^SSM lepton number is broken 
explicitly in the superpotential by terms which are trilinear as well as linear in singlet 
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neutrino superfields. In addition to that there are lepton number conserving terms involving 
the singlet neutrino superfields with dimensionless neutrino Yukawa couplings. After the 
electroweak symmetry breaking these terms can generate the effective bilinear R-parity 
violating terms as well as the AL =2 Majorana mass terms for the singlet neutrinos in the 
superpotential. In general, there are corresponding soft supersymmetry breaking terms in 
the scalar potential. Thus the parameter space of this model is much larger compared to 
the bilinear R-parity violating model. Hence, in general, one would not expect a very tight 
correlation between the neutrino mixing angles and the ratios of decay branching ratios of 
the LSP. However, under certain simplifying assumptions (as discussed in Sec. 3.2), one 
can reduce the number of free parameters and in those cases it is possible that the above 
correlations reappear. As mentioned earlier, this has been studied in great detail for the 
two body £ — W final states in [73] and for all possible two and three body final states 
in [79] . Let us note in passing that such a nice correlation is lost in the general scenario of 
bilinear-plus-trilinear R-parity violation [53]. 

Another important difference between /iz^SSM and the bilinear R-parity violating 
model in the context of the decay of the LSP (assumed to be the lightest neutralino in 
this case) is that in /ii^SSM the lightest neutralino can have a significant singlet neutrino 
(v^) contribution. In this case, the correlation between neutrino mixing angles and decay 
branching ratios of the LSP is different [73, 79] compared to the cases when the dominant 
component of the LSP is either a bino, or a higgsino or a Wino. This gives us a possibility 
of distinguishing between different R-parity violating models through the observation of 
the decay branching ratios of the LSP in collider experiments [73, 79]. In addition, the 
decay of the lightest neutralino will show displaced vertices in collider experiments and 
when the lightest neutralino is predominantly a singlet neutrino, the decay length can be 
of the order of several meters for a lightest neutralino mass in the neighbourhood of 50 
GeV [79] . This is very different from the bilinear R-parity violating model where for a Bino 
LSP of similar mass the decay length is less than or of the order of a meter or so [56] . 

The paper is organized as follows. We start with a brief introduction to the model in 
section Q and discuss the electroweak symmetry breaking conditions. The neutrino sector is 
discussed in section |3| in details, accompanied with necessary analytical results. We discuss 
the observed pattern of neutrino masses and mixing upto one-loop corrections. We present 
a comprehensive discussion on the results of our numerical analysis of neutrino masses 
and mixing in section |^. The three broad scenarios, namely, normal hierarchy, inverted 
hierarchy, and quasi-degenerate neutrinos, are taken up in turn in this section. We conclude 
in section |5|. Various technical details, such as different mass matrices, couplings, Feynman 
rules and the expressions for one-loop contributions are relegated to the appendices. 



2. Electroweak symmetry breaking in /iz/SSM 

The superpotential for ^i/SSM includes three gauge-singlet right handed neutrino super- 
fields {Of (i = e,fi,T)) along with the usual MSSM superfields. The superpotential of 
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HvSSM along the lines of ref. [71] is 



1 

\i Tja Tjb , .Ajk^-C'-C'-c 



(2.1) 



where Hi and H2 are the Higgs superfields that have Yukawa couplings with down- and 



(down-type) quark superfields, Li are doublet lepton superfields, and ij are singlet charged 
lepton superfields. The absence of any bilinear terms in the superpotential is ensured by 
imposing a Z3 symmetry (which is also used in case of NMSSM). The effective /x-term 
is given hy fi = Yl'^^'^i^ where v'l is the VEV obtained by the 'z'-th right handed sneu- 
trino (scalar component of Pf) after EWSB. The characteristic bilinear Rp violating terms 
{£iLiH2) appear in a similar way after the EWSB. These terms are given by e* = ^ Yy-' v'^. 
The last two terms of the superpotential (see eq.(|2.lD) violate Rp through L-violation. The 
last term, with the coefficient k^^^ , is included in order to avoid an unacceptable axion 
associated to the breaking of a global U(l) symmetry [83]. This term generates effective 
Majorana masses for the singlet neutrinos at the electroweak scale. 

It has been shown earlier [71-73] that the above i?p-violating superpotential (see 
eq.( |2.lD ) provides the minimal structure, sufficient for both generating a neutrino mass 
pattern, and offering a solution to the //-problem. There have been studies [84,85] on Rp- 
violating scenarios including right handed neutrinos, which use a subset of this minimal 
superpotential, but without any attempt to address the //-problem. Similar remarks apply 
to earlier works [86] aimed at explaining the baryon asymmetry of the universe through 
leptogenesis, using the term i''^HiH2- 

The frequently discussed trilinear L-violating terms, driven by the well-known A- and 
A'-type couplings, need not appear explicitly in this model. The reason is the following; 
the above superpotential, as has already been stated, can lead to 'bilinear' terms of the 
form LiH2 once the right sneutrinos acquire VEV, and, once such terms arise, they can 
effectively lead to the A and A'-type terms [30,31]. Here, as a digression, let us mention 
that the spontaneous breakdown of the symmetry through right-sneutrino VEV can in 
general lead to the formulation of domain walls [87]. The associated problems can, however, 
be ameliorated through well-known methods [88]. 

Coming back to /iz^SSM, if we confine ourselves to the framework of supergravity me- 
diated supersymmetry breaking, the Lagrangian Csoit-, containing the soft-supersymmetry- 
breaking terms is given by 



up-type quark; 



;s, respectively. Qi are doublet quark superfields, -u^ (dp are singlet up-type 





(2.2) 
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The first two lines of eq.( |2.2D consist of squared-mass terms of squarks, sleptons and Higgses, 
the next two hnes contain the trihnear scalar couplings, while in the last line, M3, M2, and 
Ml represent the Majorana masses corresponding to SU{3), SU{2) and U{1) gauginos 
A3, A2, and Ai, respectively. The tree- level scalar potential receives the usual D and F term 
contributions, in addition to the terms from £soft- 

We adhere to the CP-preserving case, so that only the real parts of the neutral scalar 
fields develop, in general, the following VEVs, 



0\ 



V2 



(2.3) 



The tree level neutral scalar potential looks like 

2 



(^^eutral) 



+ E 



i.k 



+ 



+ 



+ ( 



9i+9l 



V2 



(2.4) 



It is important to notice that the potential is bounded from below as the coefficients of 
the fourth power of all the eight superfields are positive. We further assume that all 



the parameters present in the scalar potential are real. From eq.(2^), the minimization 
conditions in terms of v\^vi^ v\ can be derived (the equations are provided in appendix 
^). The minimization conditions for //z/SSM have also been addressed in [72, 73, 79] . 
Similar conditions, but for complex VEVs, have been discussed in ref. [80]. Note that in 
order to generate correct order of magnitudes for the light neutrino masses through the 
TeV scale seesaw mechanism, one requires smaller values for neutrino Yukawa couplings 
(yj^' ~ O (10-6)) and left handed sneutrino VEVs {y\ ~ O (10^^) GeV). 



3. The neutrino sector 
3.1 Neutral fermions 

In this model, three SU{2)l doublet neutrinos (fj) and three gauge-singlet right handed 
neutrinos (u^) mix with the MSSM neutralinos (two neutral gauginos and two neutral 
higgsinos) due to L-violating interactions (see eq.( |2.lD ). The resulting neutralino mass 
matrix therefore is of dimension 10 x 10. The mixing among various current eigenstates 
are governed by the VEVs of various neutral scalar fields (namely, Hf,H2,i^i,i^f). This 
matrix has been addressed in refs. [71-73, 79] for real VEVs, and in ref. [80] for complex 
VEVs. 
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In the weak interaction basis, defined by, 

M/O^ = (5°, WlH^„ Hi vl i.^, V,, U^, Ur 
the neutral fermion mass term in the Lagrangian is of the form 



^neutral 



where Mn is the 10 x 10 modified neutraUno mass matrix, and is given by 



^7x7 "iLr 

"2-3x7 03x3 



Here, using eq.( |A.5|) 



/ Ml 




91 



M7 



x7 



91 



V2 



\ 


/i —\^V2 —\^V2 —X^V2 







p' 



-X^V2 p^ 



2uf 2uT luf 



and 



-\^V2 p^ 2uT 2uT luf 



\ -X'V2 /O^ 2n^^ 214^^^ 2uY / 



^< rl Y^^V2 Yu^V2 Y^''V2\ 



"^3x7 



^< r^' yr^2 i^rt;2 1^^"^;2 



(3.1) 



(3.2) 



(3.3) 



(3.4) 



(3.5) 



V-^< ^< Tl Y--V2 Y^^V2 Y--V2 j 

The matrix M7X7 contains the 4x4 block (upper left) of MSSM neutralinos as well as 
a 3 X 3 block (bottom right) of gauge-singlet neutrinos and mixing terms between them. 
The null 3x3 block in 7W„ signifies the absence of Majorana mass terms for the left 
handed neutrinos. The elements of m^xr contain either left handed sneutrino VEVs {v'^) 
or Higgs VEVs multiplied by neutrino Yukawa couplings iY^-'), and hence, are of much 
smaller magnitudes compared to the entries of Mjx7- This feature ensures a seesaw-like 
structure of Ain- 

This 10 X 10 symmetric matrix Ain can be diagonalized with a 10 x 10 unitary matrix 
N to obtain the physical neutralino states. The mass eigenstates are defined by, 



Xi 



Ni.^l i,j = l,...,W, 



(3.6) 
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where satisfies 

N*MnN-'=M%, (3.7) 

with the diagonal neutralino mass matrix denoted as Seven eigenvalues of this matrix 

turn out to be heavy, i.e. of the order of the electroweak scale, and thus correspond to the 
physical neutralinos. The remaining three light eigenvalues correspond to the masses of 
three SM- neutrinos. One can therefore write eq.(^.7|) alternatively as 

N*MnN-^ = diag(m~o,mj), (3.8) 

where i = 1, 7 and j = 1, 2, 3. 

Assuming small Rp violation, it is possible to carry out a perturbative diagonalization 
of the 10 X 10 neutralino mass matrix (see [89]), by defining [90] a matrix ^ as 

e = m3x7.M7-V. (3.9) 

If the elements of ^ satisfy <C 1, then this can be used as an expansion parameter to get 
an approximate analytical solution for the matrix N (see eq.(3.7)). A general expression 
for the elements of ^ with simplified assumptions can be written in the form Aai+Bbi+Cci, 
where 

a, = Yi'v2, c, = v'i, b, = {Yi'vi + 3A^;^) = (a^ cot /3 + 3Aci), (3.10) 

with i = e,fi,T = 1,2,3, tan/3 = ^ and A,B,C are complicated functions of various 
parameters of the model. The complete expressions for the elements of ^ are given in 
appendix M. Here we neglect the subdominant terms O ~ , ^^4^, where fh is the 

electroweak (or supersymmetry breaking) scale. 

The mixing matrix in leading order in ^ is given by 



n rrT V\ \ ..t • (3-11) 




The 10 X 10 neutralino mass matrix A4n can approximately be block-diagonalized to the 
form diagfMixY, M^^'^'^'^^ ), by the matrix defined in eq. (|3.11 ). The matrices Af and Ui^, 
defined in eq.( |3.11 ), are used to diagonalize AfyxT and M^^^'^'^^ in the following manner, 

AA*M7x7AAt = diag(m~o), 
UlM~U, = diag(mi,m2,m3). 

(3.12) 

3.2 Seesaw mechanism and tree level neutrino mass 

The effective light neutrino mass matrix Mj^^'^^"''^, arising via the seesaw mechanism in 
presence of explicit lepton number violation, is in general given by 

M~ = -m3x7Mf^V"^L7- (3.13) 

The eigenvalues and eigenvectors of the 3x3 matrix M^'^^^^'^ were computed in ref. [73] 
with the simplifying assumption of a flavour diagonal structure of the neutrino Yukawa 
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couplings Yl'' . Three flavour global neutrino data were fitted with this assumption, and it 
was observed that all three neutrinos acquire masses even at the tree level. 

An approximate analytical expression for the elements of Mff'^^"'^ at tree level, as 
obtained in ref. [73], with certain simplifying assumptions is given by 

{M~)ij = - 35i,) + -^b,bj. (3.14) 



One can further rewrite eg. ( 3.14 ) in an elucidate form given by 
where aj and q are given by eq.( 3.10| ) and 



with 



1 2Aifc^ 

2AXfi 2Ancot(3 
h = ^, /3 = ^^' (3-16) 



/i = ^Xv", A = {kv"'^ + XV1V2), 

A = X'^ivf + u|)2 + 4XkviV2v'''^ - 'iXA^M, 

L = A + A_ (3.17) 

M Ml Ml ^ ^ 



The reason for recasting eq.( 3.14| ) in terms of Oj and becomes clear when we will discuss 



our numerical results in section ^ In (M^'^^*"'")jj, we neglect the subdominant terms of 
the order of ^"-^^ , and Y^v' . For the convenience of the reader, let us also mention 

here that we choose A*, {AxXy, k^^^, {AkkY-'^ and all soft masses to be flavour diagonal 
and flavour blind, The neutrino Yukawa couplings (Y^-') and the corresponding soft terms 
{A^Y^y^ are, however, chosen to be flavour diagonal. 

3.3 One loop corrections to the self energies 

In the regime of renormalizable quantum field theories, stability of any tree level analysis 
must be re-examined in the light of radiative corrections. Following this prescription, 
the results of neutrino masses and mixing will be more robust, once tree level analysis 
is further improved by incorporating radiative corrections. The radiative corrections may 
have sizable effect on the neutrino data at one-loop level. Thus, although all three SM 
neutrinos acquire non-zero masses in the //z^SSM even at the tree level [73], it is interesting 
to investigate the fate of those tree level masses and mixing when exposed to one-loop 
corrections. With this in view, in this section we perform a systematic study of the neutrino 
mass and mixing with all possible one-loop corrections both analytically and numerically. 
In the subsequent sections, while showing the results of one-loop corrections, we try to 
explain the deviations (which may or may not be prominent) from the tree level analysis. 
The complete set of one- loop diagrams are shown in figur^. Before going into the details. 
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let us discuss certain relevant issues of one-loop correction and renormalization for the 
neutralino-neutrino sector. The most general one-loop contribution to the unrenormalized 
neutralino-neutrino two-point function can be expressed as 

iS|^o(p) = ^{y[PL4(/) + Pi?S,^(p2)] _ [PLU^ip^) + PnU^,{p')]}, (3.18) 

where Pl and Pr are defined in eq.( pi.2| ), i, j = 1, 10 and p is the external momentum. 
The unrenormalized self-energies S and 11 depend on the squared external momentum (p^). 
The generic self energies "E-j , H-j of the (Majorana) neutrino must be symmetric in 
its indices The resulting one-loop corrected mass matrix using dimensional reduction 
(DR) scheme [91] is given by 

(_^^ree+i-loop)., ^ ^^..(^^j^*. + ^ (uY^imj) + IjY.im^,) " m^oSj<(mf ) - m^,^im] 
with 



(3.19) 



where the tree level neutralino mass {m^o) is defined at the renormalization scale ^r, set 
at the electroweak scale. Here, the word neutralino mass stands for all the ten eigenvalues 
of the 10 X 10 neutralino mass matrix. The self-energies S, 11 are also renormalized in the 
DR scheme [91] and denoted by S and IT respectively. The detailed expressions of S]^ and 
HYj are given in appendix ^. 

3.4 Radiative corrections to neutrino mass terms 

In this section we consider the effect of radiative corrections to the light neutrino masses. 
Let us recapitulate some of the earlier work regarding one-loop corrections to the neutralino- 
neutrino sector. The complete set of radiative corrections to the neutralino mass matrix in 
the Rp conserving MSSM was discussed in ref. [92], and the leading order neutrino masses 
has been derived in ref. [6]. One- loop radiative corrections to the neutrino- neutralino mass 
matrix in the context of a -Rp-violating model were calculated in ref. [25] using 't-Hooft- 
Feynman gauge. In ref. [42], R^ gauge has been used to compute the corrections to the 
neutrino-neutralino mass matrix at one-loop level in an i?p-violating scenario. Neutrino 
mass generation at the one- loop level in other variants of i?p-violating MSSM has also been 
addressed in refs. [19,39,40,44-51,84,99]. 

We begin by outlining the strategy of our analysis. We start with a general 10 x 10 
neutralino matrix, with off-diagonal entries as well, which has a seesaw structure in the 
flavour-basis (see eq.(|3.3D). Schematically, we can rewrite eq.(|3.3|) as. 



Mn=[ ^ (3.21) 

where the orders of the block matrices are as those indicated in eq.(^.3|), and the subscript 
'/' denotes the flavour basis. Here Mf stands for the 7x7 Majorana mass matrix of 
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the heavy states, while rriDf contains the 3x7 Dirac type masses for the left handed 
neutrinos. In the next step, instead of utilising the seesaw structure of this matrix to 
generate the effective light neutrino mass matrix for the three active light neutrino species, 
we diagonalize the entire 10 x 10 matrix Ain- The diagonal 10 x 10 matrix Mfy thus contains 
tree level neutralino masses, which we symbolically write as 

M% = { ^'7 " 1 , (3.22) 




where ("T-m) are the masses of the heavy states (left handed neutrinos). At this 
stage we turn on all possible one-loop interactions, so that the 10 x 10 matrix Ai^ picks 
up radiatively generated entries, both diagonal and off-diagonal. The resulting one-loop 
corrected Lagrangian for the neutralino mass terms, following eq.( |3.2| ), can be written as 

>C' = -^X°" {Ml +M')x' + H.c, (3.23) 

where contains the effect of one-loop corrections. The 10 x 10 matrix is diagonal, 
but the matrix is a general symmetric matrix with off diagonal entries. 
One can rewrite the above equation, using eqs.(^]^ and 3.7), as 

£' = -^^^^ {Ma + N^M^N) + H.c. (3.24) 

This is nothing but the one-loop corrected neutralino mass term in the Lagrangian in the 
flavour basis. Symbolically, 

£' = -l^o^X'^o + H.c, (3.25) 



with the 10 x 10 matrix Ai' having the form 

(3.26) 



The quantities AMf and Amj stand for one-loop corrections to the heavy neutralino 
states and light neutrino states respectively, in the flavour basis The entity Atjid^. 
arises because of the off diagonal interactions, i.e. between the heavy neutralinos and the 
light neutrinos, in the same basis. Note that all of AMf, Am^)^. , Anij in the xo basis are 
given by the second term on the right hand side of eq.( |3.19| ). We suitably transform them 
into the basis with the help of neutralino mixing matrix A^. Interestingly, the matrix 
M.' once again possesses a seesaw structure, and one can therefore write down the one-loop 
corrected effective light neutrino mass matrix as 

(iVr')eff « Anif - {mof + AmDf){Mf + AMf)-\{mDj + AmOfY)- (3-27) 

Let us now present an approximate form of eq.( p.27| ). For simplicity, let us begin by 
assuming the quantities present in eq.( ^.27|) to be c- numbers (not matrices). In addition, 
assume Mf > AMf (justified later), so that eq. (|3^ ) may be written as, 
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with (5 = ^1 Mj^j • Now, even when Amof ~ rriD^. and AMf ~ Mf, eq.(3.25 

looks hke 




2567r4 V M 



ruD, 



2 



(3.29) 



Thus, up to a very good approximation one can rewrite eq. (|3.29D as 

{M'^),s^Amf-Mf^-j^j . (3.30) 
Reimposing the matrix structure and using eq.(3.13), eq.(3.30) can be modified as, 



{M" )eff ^ Amj + M""". (3.31) 

The eigenvalues of the 3x3 one-loop corrected neutrino mass matrix {M^ )eff thus cor- 
respond to one-loop corrected light neutrino masses. In conclusion, it is legitimate to 
calculate one-loop corrections to the 3x3 light neutrino mass matrix only, and diagonalize 
it to get the corresponding one- loop corrected mass eigenvalues. 

Let us denote the one-loop corrections to the masses of heavy neutralinos and light 
neutrinos in the basis by AM and Am respectively. The one-loop corrections arising 
from neutralino-neutrino interactions is denoted by AniD in the same basis. The tree level 
neutralino mixing matrix N can then be written as, 

N=( ^7x7iV7x3y (332) 

where the entries of the matrices N^xs, ^3x7 are ~ O (m'^/M^o), due to very small 
neutrino-neutralino mixing [93]. The quantities and M^o represent the Dirac mass of 
neutrino and the Majorana mass of neutralino. From eq.( p.24| ), it is easy to figure out the 
relation between Am and Amj as, 

Amj = iVf^3AMiV7x3 + iVfxsAml^iVgxs + Ng-^^AmoNry^s + Ng-^^AmN^y.^. (3.33) 

Typically, for a Dirac neutrino, the mass is ^ O (10~^ GeV), while for a neutralino, the 
mass is ~ O (10^ GeV). This means that the entries of the off-diagonal blocks in eq.( 3.32| ) 



are ^ O (10 ^). Therefore, for all practical purposes, one can neglect the first three terms 
in comparison to the fourth term on the right hand side of eq.( |3.33| ). Thus, 

Amf^ N^x^AmN^xs- (3.34) 

up to a very good approximation. With this in view, our strategy is to compute the one- 
loop corrections in the basis first, and then use eq.( |3.34 ) to obtain the corresponding 



corrections in the flavour basis. Finally, we diagonalize eq. (|3.31| ) to obtain the one-loop 
corrected neutrino masses. We have performed all calculations in the 't-Hooft-Feynman 
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gauge. Let us also note in passing that the form of eq.( |3.19| ) predicts off-diagonal entries 
{i 7^ j). The off-diagonal elements are responsible for the mixing between diagonal entries, 
they become dominant only when (jn^o — m^o^ ^ ^ some electroweak scale mass, 

and then, one can choose p'^ = m? = {m?^i3 + m?~o)/2 for external momentum [251. Thus, 
one can conclude that unless the tree level masses are degenerate, the off-diagonal radiative 
corrections can be neglected for all practical purposes, when at least one index i or j refers 
to a heavy states. 

The self-energy corrections contain entries of the neutralino mixing matrix through 
the couplings O^^''' (see, appendix E). This is because, the self energies Sy and i7jj in 
general contain products of couplings of the form Oj . The matrix N, on the 



other hand, contains the expansion parameter ^ in the leading order (see eg. ( 3. 111) ). This 



observation, together with the help of eq. ([B.l|) , help us to express the effective structure 



of the one-loop corrected neutrino mass matrix as, 

[{M''')eff]ij = Aimaj + A2CiCj + A'i{aiCj + OjCj), (3.35) 

where and are given by eq.( |3.10| ) and AiS are functions of our model parameters and 
the Passarino-Veltman functions {Bq^Bi) [94-96] defined in appendix The form of the 
loop corrected mass matrix thus obtained is identical to the tree level one (see, eq. (|3.15|) ) 
with different coefficients Ai, A2 and A-^ arising now. 

Note that the one-loop diagrams in figure||, contributing to the neutrino mass matrix 
are very similar to those obtained in bilinear R-parity violating scenario [42-45,50,51]. 
However, it has been pointed out in Ref. [79], that there is a new significant contribution 
coming from the loops containing the neutral scalar and pseudoscalar with dominant singlet 
component. This contribution is proportional to the mass-splitting between the singlet 
scalar and pseudoscalar states [97-99]. The corresponding mass splittings for the doublet 
sneutrinos are much smaller [79]. 

3.5 Neutrino mixing 

The unitary matrix which diagonalizes the 3x3 light neutrino mass matrix, can be 
parametrized as [100], 



C12C13 



S12C13 



Sl3 



-S12C23 - C12S23S13 C12C23 - S12S23S13 S23C13 



\ S12S23 - C12C23S13 -C12S23 - S12C23S13 C23C13 / 



(3.36) 



provided the charged lepton sector is in the mass-basis (see appendix ^). Here c 



cos ( 



sinOij, and all CP violating phases (Dirac or Majorana) are set to zero. 
The experimental data on neutrino oscillations [101, 102] indicate that a bilarge pattern 
of mixing of neutrinos may be preferred, which means that two of the mixing angles must 
be large. As a first approximation, one can work with ^23 = 45°, sin0i2 = and 
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Figure 1: One-loop diagrams contributing to the neutralino masses. The various contributions are 
arising from (clockwise from top left) (a) neutralino-neutralino-neutral scalar loop, (b) neutralino- 
neutralino-neutral pseudoscalar loop, (c) neutralino-neutralino-Z^ loop, (d) neutralino-chargino- 
charged scalar loop, (e) neutralino-chargino-W^ loop, (f) neutralino-quark-squark loop. 

0°, which is often referred to as the 'tribimaximal structure' [103]. Following the 
above discussions, we can write down the tree level PMNS matrix as, 

= ivf,3, (3.37) 

where A'^xs matrix that diagonalizes the tree level mass matrix M^^'^^"''^ (eq.( |3.1t^ )). 

Similarly, the 3x3 unitary matrix that diagonalizes the one-loop corrected neutrino 
mass matrix (M^ (eq.(|^)), can be denoted as Ul^. Symbolically 

Z^r'(^"')cff^^^ = diag(m;, m'^, m'g), (3.38) 

with m'l, 771-2, "^3 ^ three one- loop corrected light neutrino masses. A similar relation 
for the tree level calculation is given by the second equation of ( p.l2| ), with m^'s as the tree 
level neutrino masses. 

When we include one-loop corrections, the PMNS matrix Ul is defined as 

K = A^x3, (3.39) 

where N^^,^ is the matrix which diagonalizes the one- loop corrected mass matrix (M^')eff. 
The scheme for obtaining the loop-corrected neutrino mixing matrix is shown in figure^. 

4. Numerical results of neutrino mass and mixing 

In this section we present and explain the results of our numerical analysis. Let us begin 
with a brief outline of the different well known schemes of light neutrino masses as favoured 
by experiments. These are (i) The normal hierarchy: mi < m2 ~ ^ ^"^Liar ' "^3 ~ 
^y\Aml^^\, (ii) The inverted hierarchy: mi m2 ~ y^\Arn^^\, nis ^ y/\Am^^\ and (iii) 
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Flow chart of our calculational approach 



Figure 2: Flow chart of the calculational procedure that we have followed in order to calculate 
the neutrino masses and mixing at the tree level and at one-loop. Here U°pmn°s'^ ~ ^'v ^pmns — 

7 J ntree rmass T\/rtree l\/f seesaw none— loop rl i\/rone — loop fl\/fJ^'\ 

"i'' '-flavour-basis — '-neutral^ ^^'seesaw — ^"u t — >- ^^'seesaw — Jeff m 

the text. 



The quasi-degenerate pattern: mi ~ m2 ~ S> \/\Am'^^^\, where mi, m2 and are 
the three light neutrino masses. Here, Am^^^^^ = — and Am^^^ = — m^mf) for 
Normal (inverted) Hierarchy. The possibility of more than one scheme of neutrino masses 
essentially stems from our lack of knowledge of the signs of the squared mass differences, 
or the value of individual masses. Such being the case, the numerical analysis must be 
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subjected to address all the three probable schemes. 

The numerical calculations have been performed with the help of a code developed 
by us in Fortran and the results have been cross-checked using another Mathematica [104] 
based program we developed. In the code, we keep the left handed sneutrino VEVs {v'), 
the right handed sneutrino VEVs (v^) and the neutrino Yukawa couplings (Y^) as free 
input variables, and then scan over the parameter space for a region consistent with the 
three flavour global neutrino data. For all our numerical analysis we keep the right handed 
sneutrino VEVs fixed at some chosen values (see table |l]). The only exception to this is 
when we study the correlation of neutrino data with the bilinear iZp-violating parameter 
Ei = Yi^-'vj. For that particular study, we vary the right handed sneutrino VEVs (consistent 
with EWSB conditions) within the mass scale -895 GeV to -565 GeV. The table I shows 
our choice of the sample parameters for the numerical analysis. The relation between the 
gaugino soft masses Mi and M2 are assumed to be GUT (grand unified theory) motivated, 
so that, at the electroweak scale, we have Mi : M2 = 1:2. We choose Mi = 110 GeV. 



tan/3 


A 


AxX 








{m%) 




10 


0.10 


-ITeV X A 


0.45 


ITeV X K 


ITeV X Yjj 


3002 QgY2 


-595 GeV 



Table 1: Choice of parameters for numerical analysis consistent with the EWSB conditions. 



We scanned the parameter space comprising of the left handed sneutrino VEVs and 
the neutrino Yukawa couplings extensively and found certain ranges of these parameters 
appropriate for various hierarchical schemes of the light neutrino masses. In table g some 
sample values for these six parameters for the different mass schemes of neutrinos are 
given. These values are just for illustration and the parameters were scanned around these 
numbers to generate the plots shown in this section. The other relevant parameters have 
values as mentioned in table |l[ 

While fitting the three flavour global neutrino data, we consider constraints arising 
from the oscillation data as well as from the non-oscillation data. We probe the effects of 
these constraints for both tree level and (tree -|- one-loop) level analyses. The oscillation 
data constrain the solar and atmospheric mass squared differences, namely, Am^^^^j, and 
Am^j^, and three neutrino mixing angles ^13, 9i2, ^23- The present 3a limits are [2,3], 



7.05 X lO-^eV^ < Amli,, 
2.07xlO-3eV2<|Am2,„| 



< 8.34 X 10~^eV^ 

< 2.75 X lO^^eV^ 



0.25 < sin^ 6^12 < 0.37, 
0.36 < sin^ 6^23 < 0.67, 
sin^ 013 < 0.056. 



(4.1) 



The non-oscillation constraints follow from experiments like f3 decay [105-107], neutrinoless 
double beta decay (Oz^/3/3) [108-110] (this is also sensitive to Majorana nature and phases), 
and from cosmology [111]. Here we set all the Majorana phases to be zero, as we are 
dealing with a CP— preserving situation. 
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X lO^(GeV) 


V 


V 










Normal hierarchy 


3.550 


5.400 


1.650 


0.730 


10.100 


12.450 


Inverted hierarchy 


12.800 


3.300 


4.450 


8.350 


8.680 


6.400 



Table 2: Values of neutrino Yukawa couplings and left handed sneutrino VEVs, used as sample 
parameter points for numerical calculations. These are the values around which the corresponding 
parameters were varied. 



The set of non-oscillation constraints are given as, 



= <i.80eV (/3 decay), 

0.00 eV < rrififi = \ ^Uujmi\ < 0.25 eV (Oz^/3/3), 

^^mj < 1.30 eV (cosmology), 



(4.2) 



where mjS are three light neutrino masses, and U^^^s are the elements of the first row of 
the neutrino mixing matrix (see eq.( |3.36| )). 

4.1 Normal hierarchy 

In the normal hierarchical pattern of the three light neutrino masses, the atmospheric and 



2 

atm 



-ml and Am^^,^^ 



nil, 



the solar mass squared differences, given by Am 
largely governed by the higher mass squared in each case, namely, m| and m|, respectively. 
Before going into the discussion of the variation of the mass-squared values with the model 
parameter, some general remarks are in order. First of all, note that in eq.( 3.1C| ), if we 

Cj [80]. Thus, for large f we remind ourselves that 



choose such that ^ ' then h. 



in eq.( |3.14 ), the effective light neutrino mass matrix has two types of seesaw structures [73]. 
The first one is the ordinary seesaw, given by 

~ (4.3) 

where Oj = Y^'^V2 represents the Dirac mass term for neutrinos, and m^^c = 2kv^ stands 
for the Majorana mass term of the right handed neutrino. The second type is called the 
gaugino seesaw, in which the role of the Dirac mass terms are played by giQ and g2Ci, 
where gi, g2 are the U{1) and the SU(2) gauge couplings respectively and q stands for the 
left handed sneutrino VEV v'^. The role of the Majorana masses are played by the gaugino 
soft masses Mi, M2. This seesaw relation is given as 



Ml 
M' 



M2 



(4.4) 



where the subscript 



1,2,3 = e, fi, T and M is the reduced gaugino mass defined by 



M here plays the role of the effective heavy mass provided by the neutral electroweak 
gaugino sector, and the effect is closely analogous to Type-Ill seesaw mechanism [112]^. 
As discussed after eq.( p.35 ), when one-loop corrections are added, the neutrino masses are 
still determined by the quantities af and c? bf for large v'^. 
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Figure 3: Neutrino mass squared values (mf) vs (left panel) and vs (right panel) plots 
for the normal hierarchical pattern of light neutrino masses, i = e, ^, t. 

In the subsequent plots, we show the variation of the neutrino squared masses (m?) 
and the atmospheric and solar mass squared differences with the square of the seesaw 
parameters and . Results are shown for the tree level as well as the one- loop corrected 
neutrino masses. These plots also demonstrate the importance of one-loop corrections to 
neutrino masses compared to the tree level results. 

""^We thank Anjan Joshipura for pointing this out to one of the authors in a private discussion. 
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Typical mass spectra are shown in figure ^. Note that a particular model parameter 
has been varied while the others are fixed at values given in tables |l| and ^. The effective 
light neutrino mass matrix given in eq.(3.14) suggests that as long as v[ ^> ~^x^ ^"^^ 
K S> A, the second term on the right hand side of eq.( |3.14 ) dominates over the first term 
and as a result the heaviest neutrino mass scale {m^) is controlled mainly by the gaugino 
seesaw effect. This is because in this limit hi ~ Cj, and, as discussed earlier, a neutrino 
mass matrix with a structure (TOv)ij ~ can produce only one non-zero neutrino mass. 
This feature is evident in figurejs], where we see that m| increases as a function of cf/M^. 
The other two masses are almost insensitive to cf/M. A mild variation to comes from 
the combined effect of gaugino and ordinary seesaw. On the other hand, the two lighter 
neutrino mass scales (m^ and mf) are controlled predominantly by the ordinary seesaw 
parameters af/rriyc. This behaviour is observed in the right panel figures of figure^. The 
heaviest neutrino mass scale is not much affected by the quantities af/niyc. 

One can also see from these plots that the inclusion of one-loop corrections, for the 
chosen values of the soft SUSY breaking parameters, reduces the values of and mf, 
while increasing the value of only mildly. This is because, with such a choice, the one- 
loop corrections cause partial cancellation in the generation of mi and m2- For the heaviest 
state, it is just the opposite, since the diagonalization of the tree-level mass matrix already 
yields a negative mass eigenvalue, on which the loop correction has an additive effect. If, 
with all other parameters fixed, the signs of A and A\ are reversed (leading to a positive 
// in the place of a negative one), mi, ra2 and are all found to decrease through loop 
corrections. A fiip in the sign of k and the corresponding soft breaking terms, on the other 
hand, causes a rise in all the mass eigenvalues, notably for mi and m2. 

In the light of the discussion above, we now turn to explain the variation of Am^^^ and 
^fn^soiar with cf/M"^ and af/mlc shown in figure^and figure|5|. For our numerical analysis, 
in order to set the scale of the normal hierarchical spectrum, we choose m2\max < 0.011 eV. 
The left panel in figure^ shows that Am^^^ increases more rapidly with c^^/M^, whereas 
the variation with d^/M'^ is much slower as expected from figure^ Similar behaviour is 
shown for the one- loop corrected Am^^^. The small increase in the one- loop corrected 
result compared to the tree level one is essentially due to the splitting in m| value as 
shown earlier. The variation of Am^^j^j. with c|/M^ can be explained in a similar manner. 
Obviously, in this case the one-loop corrected result is smaller compared to the tree level 
one (see, figure^). However, one should note that Am^^^^^ falls off with c^/M^ as opposed 
to the variation with respect to the other two gaugino seesaw parameters. This is due 
to the fact that slightly decreases with c^/M^ but show a slow increase with respect 
to Cg/M^ and c^/M^. The dark solid lines in all these figures show the allowed values of 
various parameters where all the neutrino mass and mixing constraints are satisfied. 

The variation of Am^j^ and Am^^^^^ with af/rri^c in figure^ can be understood in 
a similar way by looking at the right panel plots of figur^. Am^^^ shows a very little 
increase with flg^^/m^c as expected, whereas the change is more rapid with a^/m^c for the 
range of values considered along the x-axis. As in the case of figure^ the solid dark lines 
correspond to the allowed values of parameters where all the neutrino mass and mixing 
constraints are satisfied. 
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Figure 4: Atmospheric and solar mass squared differences (Ato^j,„, Am^^j^j^) vs plots for the 
normal hierarchical pattern of light neutrino masses, i = e, /i, r. The full lines are shown for which 
only the constraints on Ato^^^j^j, is not within the 3cr limit. The dark coloured portions on these 
lines are the values of parameters for which all the neutrino constraints are within the 3(t limit. The 
red (yellow) coloured lines in the plots correspond to the tree (one-loop corrected) regions where all 
the constraints except Am^^j^j, are within Scr allowed region. Other parameter choices are discussed 
in the text. 



For higher values of a^^/m^c, m\ increases very slowly with these parameters (see, 
figure^) and this is reflected in the right panel plots of figure^, where Am^^^^^ shows a very 
slow variation with Og ,-/m^c. On the other hand, increases more rapidly with a^/m^c, 
giving rise to a faster variation of ^rn^goiar ■ plots of figur^ show that larger values of 
Yukawa couplings are required in order to satisfy the global three flavour neutrino data, 
when one considers one-loop corrected neutrino mass matrix. However, there are allowed 



-20- 



0.0026 



*> 0.0025 

0) 



I 0.0024 



NORMAL HIERARCHY 



o.oo: 



TREE + ONE LOOP 

\ 




I TREE LEVEL (all constraints within 3g limit) 
I TREE + ONE LOOP (all constraints within 3o limit) 



9e-05r 



> 8e-05 



E 7e-05 

< 



NORMAL HIERARCHY 



05 4e-05 5e-05 6e-05 

4,2, ,,2. 

a / m c (eV ) 




TREE + ONE LOOP 
\ 



. ■ TREE LEVEL (all constraints within 3a limit) 
„ _ _ m TREE 4- QNE LOOP (all constraints witl^in 3o limit) 

7e-05 °® 3i-05 4e-05 5e-05 6e-05 7e-05 

4,2, ,,2. 

a / m c (eV ) 



0.0026 r 



"> 0.0025 

■Si. 

E 

CM 

I 0.0024 



NORMAL HIERARCHY 



TREE + ONE LOOP 




TREE LEVEL (all constraints within 3a limit) 
TREE + ONE ^.OOP (all constraintp within 3o limit) 



9e-05 



> 8e-05 

0) 



E 7e-05 

< 



0.0002 0.0003 

a" / c (eV^) 



0.0004 



6e 



NORMAL HIERARCHY 



■ TREE LEVEL 
(all constraints within 3a limit) 

■ TREE + ONE LOOP 
(all constraints within 3a limit) 




(9.Sooi 



0.0002 0.0003 

a" / c (eV^) 



0.0004 



0.0026 



"> 0.0025 



I 0.0024 



NORMAL HIERARCHY 



0.0023, 



TREE + ONE LOOP 




9e-05r 



> 8e-05 



E 7e-05 



NORMAL HIERARCHY 







I TREE LEVEL (all constraints within 3a limit) 
I TREE + ONE LOOP (all constraints within 3a limit) 

1e-06 2e-06 3e-06 4e-06 5e-06 

a /m . (eV ) 



6e-05, 



'0 



TREE + ONE LOOP 



I TREE LEVEL (all constraints within 3a limit) 

I TREE + ONE LOOP (all constraints within 3a limit) ' 

1e-06 2e-06 3e-06 4e-06 5e-06 



4 , 2 

a / m I 



(eV" 



Figure 5: Atmospheric and solar mass squared differences (Am^j„, Am^^;^,.) vs a|/m^c plots for 
the normal hierarchical pattern of light neutrino masses with i = fi, r. Colour specification is 
same as described in the context of figure ^ Other parameter choices are discussed in the text. 

ranges of the parameters af/m^c, where the neutrino data can be satisfied with both tree 
and one- loop corrected analysis. 

We have also considered the variation of light neutrino mass squared differences with 
the effective bilinear Rp violating parameter, Si = Y^^Vj. For this particular numerical 
study we vary both Y^' and the right handed sneutrino VEVs simultaneously, in the 
suitable ranges around the values given in table |^ and Am^^^ is found to increase with 
Ej, whereas the solar mass squared difference decreases with increasing Sj. The 3a allowed 
region for the solar and atmospheric mass squared differences were obtained for the lower 
values of e^s. In addition, we have noticed that the correlations of Am^^^ with ej is sharper 
compared to the correlations seen in the case of Am^^;^^. 
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Next let us discuss the dependence of Am^j„ and Am^^^^^ on two specific model 
parameters, A and k, consistent with EWSB conditions. The loop corrections shift the 
allowed ranges of n to lower values with some amount of overlap with the tree level result. 
On the other hand, the allowed ranges of A shrinks towards higher values when one-loop 
corrections are included. These results are shown in figure^. We note in passing that the 
mass of the lightest CP-even scalar decreases with increasing A. For example, A = 0.15 can 
produce a lightest scalar mass of 40 GeV, for suitable choices of other parameters. This 
happens because with increasing A, the lightest scalar state picks up more and more right 
handed sneutrino admixture. 
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Figure 6: Plots showing the variations of Ato^j„j, Am^^;^,, with model parameters A and n for 
normal hierarchy. Values of all other parameters are given by table |l| and |^. Colour specification 
is same as described in the context of figure ^. Other parameter choices are discussed in the text. 

Finally, we will discuss the tan/3 dependence of Am^^^ and Am^^^^j,. These plots 
are shown in figure 0. The quantity Am^^^ decreases with the increasing values of tan/3 
and nearly saturates for larger values of tan/3. However, the one- loop corrected result for 
Am'^ffyi is not much different from that at the tree level for a particular value of tan /3. 
On the other hand, the solar mass squared difference initially increases with tan/3 and 
for higher values of tan P the variation slows down and tends to saturate. The one-loop 
corrections result in lower values of Am^^^^j^, for a particular tan/3. The darker and bigger 
points on both the plots of figure]^ are the allowed values of tan /3, where all the neutrino 
experimental data are satisfied. Note that only a very small range of tan/3 (~ 10-14) is 
allowed. This is a very important observation of this analysis. 

Next we will discuss the light neutrino mixing and the effect of one-loop corrections 
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Figure 7: Am^j^, Am^^;^^ vs tan/3 plots for the normal hierarchical pattern of light neutrino 



masses. The allowed values of tan/3 are shown by bold points. Other parameter choices are shown 
m table |l| and |[ 



on the mixing angles. It was shown in ref. [73] that for the normal hierarchical pattern of 
neutrino masses, when the parameter bi ~ a^, the neutrino mixing angles 623 and ^13 can 
be written as (with the tree level analysis), 

and 

On the other hand, the mixing angle 612 is a much more complicated function of the 
parameters bi and Oj and we do not show it here. Now, when bi ~ aj, we can easily see 
from eq.( |3.1[l| ), that 

t;^~^(tan/3-l). (4.8) 
This implies that for tan/3 ^ 1 (recall that the allowed range of tan/3 is ~ 10-14), 

As we have discussed earlier, for such values of v'^, the quantities 6j w Cj. Hence, the mixing 
angles 623 and ^13 can be approximately written as 

sin2^23~^^, (4.10) 



and 



,2 



sm'ei3-^^. (4.11) 
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Naively, one would also expect that sin^ 0i2 should show some correlation with the quantity 
cl/cjj^. However, as mentioned earlier, this is a very simple minded expectation since sin^ 612 
has a more complicated dependence on the model parameters. 

The variation of all three mixing angles with the corresponding parameters are shown 
in figur^. Note that in order to generate these plots, we vary only the quantities Cj and all 
the other parameters are fixed at the values given in tables |l] and ^. We have chosen the 
range of parameters in such a way that the 3-flavour global neutrino data are satisfied. The 
mixing angles have been calculated numerically by diagonalizing the neutrino mass matrix 
m eq. ( PD and in eq. (|3.35 |). As expected from our approximate analytical expressions, 
these plots show very nice correlations of the mixing angles 623 and ^13 with the relevant 
parameters as discussed in eqs.( 4.10 ) and ( [1.11 ). For example, note that when ^ Cr, 
sin^ 023 is predicted to be ^ 0.5 and that is what we observe in the tree level plot in figure^. 
However, when one-loop corrections are considered, the value of sin^ 623 is predicted to be 
somewhat on the lower side of the 3a allowed region. This can be understood by looking at 
the left panel plots of figure^, where one can see that the one-loop corrected results prefer 
lower values of and higher values of c^. Obviously, this gives smaller sin^023- On the 



other hand, the tree level analysis prefers higher values of and both lower and higher 
values of c^. This gives rise to large as well as small values of sin^^ 023- 
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Figure 8: Variation of sin^6'23 with 



sin^0i2 with sin^^ia with 



hierarchy of light neutrino masses. Other parameter choices are discussed in the text. 



for normal 



If one looks at the plot of sin^ ^13 in figure then it is evident that the amount 
fe flavour in the heaviest state (1^3) decreases a little bit with the inclusion of one-loop 

2 

corrections for a fixed value of the quantity ^^j*^^^^ . Very small sin^ ^13 demands Cg <C 



c^, c^. This feature is also consistent with the plots in figure^ The correlation of sm"^ 6 



12 
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with the ratio c^/cj^ is not very sharp as expected from the discussion given above. However, 
a large 9i2 mixing angle requires a larger value of this ratio. The effect of one-loop correction 
is more pronounced in this case and predicts a smaller value of sin^ 9i2 compared to the 
tree level result. There is no specific correlation of the mixing angles with the quantities 
af and we do not show them here. 



4.2 Inverted hierarchy 

In this subsection we perform a similar numerical analysis for the inverted hierarchical 
scheme of three light neutrino masses. Recall that for the inverted hierarchical pattern 
of light neutrino masses, the absolute values of the mass eigenvalues are such that 1712 > 
nil ^ "1-3- Thus the solar and the atmospheric mass squared differences are defined as 
Am'^f^ = ml — mg and Am^^^^^ = — m^. In order to generate such a mass pattern, the 
choices of neutrino Yukawa couplings yj* and the left-handed sneutrino VEVs v[ are shown 
in table However, these are just sample choices and other choices also exist as we will 
see during the course of this discussion. The choices of other parameters are shown in table 
|l[ The effect of one-loop corrections to the mass eigenvalues are such that the absolute 
values of masses ms and mi become smaller whereas 7712 grows in magnitude. This effect of 
increasing the absolute value of 1112 while decreasing that of rrii makes it extremely difficult 
to account for the present 3(7 limits on Am^^^^^. 

Typical mass spectra are shown in figure ^. Once again note that a particular model 
parameter has been varied while the others are fixed at values given in tables |l] and ^. As it 
is evident from these plots, the masses mi and ms are controlled mainly by the parameters 
af /m'^ , whereas the mass m2 is controlled by the seesaw parameters cf /M though there is 
a small contribution coming from af/mf, as well. 

Let us now turn our attention to the variation of |Am^^„| and Aml^i^j. with cf/M"^ and 



al jraiic shown in figure |10| and figure For our numerical analysis, we have set the scale 



of m3 as Imgl < 0.011 eV. The left panel in figure |l^ shows that |Am^j^| increases with 



c^^^jM and decreases with c^jM . This is essentially the behaviour shown by mi with 
the variation of cI/M"^. Similar behaviour is obtained for the one- loop corrected Am^^^. 
The decrease in the one-loop corrected result compared to the tree level one is due to the 
splitting in m\ value as shown in figure|9[ 

The variation of Am^^^^^ with cf/M'^ can be understood in a similar manner by looking 
at figure|9[ As explained earlier, in the case of Am^^^^^, the one-loop corrected result is 
larger compared to the tree level one. The range of parameters satisfying all the three 
flavour global neutrino data are shown by the fewer dark points on the plots. Note that 
the increase of Am^^^^^ at the one-loop level is such that we do not even see any allowed 
range of parameters when looking at the variation with respect to Cg^/M^. Once again, 
the behaviour of |Am^^^| and Am^^^^^ with the change in the parameters af/m'^c (shown 



in figure 11) can be explained by looking at the right panel plots of figure 



We have also investigated the nature of variation of | Am^^^ | and Am^ , with ef , the 



''atml O'J^i^ " solar '^ii 

squared effective bilinear i?p-violating parameters. |Am^^^| was found to increase with ef 



(the increase is sharper for e^), whereas Am^^^^,, initially increases very sharply with 
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Figure 9: Neutrino mass squared values {mf) vs (left panel) and vs ^^i- (right panel) plots 
for the inverted hierarchical pattern of light neutrino masses, i = e,/x,T. 



(particularly for and e^) ^-nd then becomes flat. In the one- loop corrected results we do 
not find any range of values for parameters where the neutrino data are satisfied. 

The variation of mass squared differences with A and k have also been analyzed. The 
variation of | Am^^„ | and Am^^^^^ with A and k are found to be opposite to those of normal 
hierarchical scenario. The one-loop corrected results do not show any allowed ranges of A 
and K (for the chosen values of other parameters) where the neutrino data can be satisfied. 

The tan j3 dependence of | Am^j„| and Am^^^^^ is shown in figurel2, One can see from 
these two figures that |Am^j„| initially increases and then starts decreasing at a value of 
tan/3 around 10. On the other hand, Am^^^^^ initially decreases and then starts increasing 
around the same value of tan/3. Note that the one-loop corrected result for |Am^j^| is 
lower than the corresponding tree level result for tan /3 < 10 whereas the one-loop corrected 
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Figure 10: Atmospheric and solar mass squared differences (|Ato^(^|, Ato^j^j^,,) vs plots for 
the inverted hierarchical pattern of light neutrino masses with i = e, /i, r. Colour specification is 
same as described in the context of figure 0. Choices of other parameters are discussed in the text. 



result for Am^^^^^ is lower than the corresponding tree level result for tan /3 > 10. For the 
chosen values of other parameters we see that the one-loop corrected analysis does not 
provide any value of tan /3 where the neutrino data can be satisfied. 

We conclude the discussion on inverted hierarchy by addressing the dependence of 
neutrino mixing angles with the relevant parameters. In figure 13 we show the variation 
of the neutrino mixing angles with the same set of parameters as chosen for the normal 
hierarchical scenario. We notice that for inverted hierarchy the quantity sin^ 023 decreases 
with increasing , which is just opposite to that of the normal hierarchy (see, figuraS^. 

1 

Nevertheless, the correlation of sin Q^j, with , ^ s is as sharp as in the case of normal 
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Figure 11: Atmospheric and solar mass squared differences (|Ato^j„|, Ato^^j^^) vs aj/m1a plots 
for the inverted hierarchical pattern of light neutrino masses with z = e, /i, t. Colour specification is 
same as described in the context of figure 0. Choices of other parameters are discussed in the text. 



hierarchy. A similar feature is obtained for the variation with j^^r^^- 

On the other hand, the correlations of sin^ 9i2 with % and % and the correlations 

2 2 '^M 

of sin^^icj with , and , are not very sharp and we do not show them here. 

There are allowed values of relevant parameters where all neutrino data can be satisfied. 
Remember that, for the plots with CjS, we varied all the CjS simultaneously, keeping the 
values of a^s fixed at the ones determined by the parameters in table |^. Similarly, for the 
variation of a^s, the quantities CjS were fixed. The inclusion of one-loop corrections restrict 
the allowed values of parameter points significantly compared to the tree level results. 
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Figure 12: \AmltJ 



^^ioiar tan/3 plots for the inverted hierarchical pattern of light neutrino 



masses. Choices of other parameters are discussed in the text. 



0.7r 



INVERTED HIERARCHY 



0.7r 



INVERTED HIERARCHY 



0.6 



0.5 

(55 



0.4 



I TREE LEVEL (aii constraints witiiin 3a limit) 

I TREE + ONE LOOP (all constraints within 3a limit) 



0.6 



'^c 0.5 
in 

0.4 



3.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

C / (c + c ) 



0.3L 



I TREE LEVEL (aii constraints within 3a iimit) 

I TREE + ONE LOOP (aii constraints within 3a limit) 



0.1 0.2 0.3 0.4 0.5 0.6 

2 , / 2 2 , 
a / (a + a ) 



0.7 



Figure 13: Variation of sin 6*23 with and ^^.^^ for inverted hierarchy of light neutrino 

masses. Choices of other parameters are discussed in the text. 



4.3 Quasi-degenerate spectra 

The discussion on the light neutrino mass spectrum remains incomplete without a note 
on the so-called "quasi-degenerate" scenario. A truly degenerate scenario of three light 
neutrino masses is, however, inconsistent with the oscillation data (see eq.( [4.1[) ). Hence, 
the quasi-degenerate scenario of light neutrino masses is defined in such a way that in this 
case all the three individual neutrino masses are much larger compared to the atmospheric 
neutrino mass scale. Mathematically, one writes mi ~ m2 ~ ms ^ yfAm^^^. Obviously, 
the oscillation data suggest that even in such a situation there must be a mild hierarchy 
among the degenerate neutrinos. 

In this section we have shown that the huge parameter space of /iz^SSM always leaves us 
with enough room to accommodate quasi-degenerate spectrum. For our numerical analysis, 
we called a set of light neutrino masses to be quasi-degenerate if the lightest among them 
is greater than 0.1 eV. We choose two sets of sample parameter points which are shown in 
table ^ (values of other parameters are same as in table ||) . For these two sets of neutrino 
Yukawa couplings (1^*) and the left-handed sneutrino VEVs {v'^) we observe the following 
patterns of light neutrino masses at the tree level 
(i) Quasi-degenerate-I: ma ^ m2 ^ mi » y^\Amlf^\, 
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yj* X lO'^ 


v\ X lO^(GeV) 


V 


V 


1/ 


A 






Quasi-degenerate-I 


19.60 


19.94 


19.99 


9.75 


10.60 


11.83 


Quasi-degenerate-II 


18.50 


18.00 


18.00 


9.85 


10.50 


10.10 



Table 3: Values of neutrino Yukawa couplings and left handed sneutrino VEVs, used as sample 
parameter points to obtain quasi-degenerate light neutrino spectrum. Around these values, the 
corresponding parameters were varied for the plots shown in figure p^. 
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(ii) Quasi-degenerate-II: m2 ^ m\ ^ ms » Am 
For case (i), we have varied the parameters around the values in table |3| and identified a 
few extremely fine-tuned points in the parameter space where either the tree level or the 
one-loop corrected result is consistent with the three flavour global neutrino data. Two 
representative spectrum as function of and are shown in figure 14. 
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Figure 14: Neutrino mass squared values (mf 



M2 



(left panel) and vs (right panel) plots 



for the quasi-degenerate pattern of light neutrino masses. 

As mentioned earlier, one can play with the model parameters and obtain a spectrum 
with a different ordering of masses termed as "Quasi-degenerate-II" in table ^ However, 
for such an ordering of masses, we found that it was rather impossible to find any region of 
parameter space where the one-loop corrected result satisfies all the constraints on neutrino 
masses and mixing. Nevertheless, we must emphasize here that it is not a completely 
generic conclusion and for other choices of soft SUSY breaking and other parameters it 
could be possible to have a spectrum like that shown in "Quasi degenerate 11" with neutrino 
constraints satisfied even at the one-loop level. On the other hand, there exists regions 
where neutrino data are satisfied at the tree level with this ordering of masses. 



5. Summary and Conclusion 

In this paper, we have performed a systematic study of neutrino masses and mixing in a 
non-minimal extension of the MSSM, known as /iz^SSM, with the complete set of one-loop 
radiative corrections, over and above the effects from the seesaw mechanism (of Types I and 
III). A set of right chiral neutrino superfields, introduced in this model provide a solution to 
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the //-problem. Lepton number is broken by one AL = 1 term comprising the right-chiral 
neutrinos and the Higgs superfields and another AL = 3 term involving only the right 
chiral neutrinos in the superpotential. The corresponding soft SUSY breaking terms in 
the scalar potential together with the F-term contribution from the superpotential, induce 
VEVs of sneutrinos of both right and left chirality. The right chiral neutrinos, together 
with neutralinos, are instrumental in the generation of light neutrino masses at the tree 
level. 

In [73], where a tree-level analysis of neutrino masses and mixing was carried out, 
it was shown that tree level masses could be generated in this fashion for all three light 
neutrinos, even with a flavour diagonal structure of neutrino Yukawa couplings {Y^-'). In 
this work we have improved the analysis with the inclusion of one-loop radiative corrections, 
for various patterns of light neutrino masses, namely, the normal, inverted and quasi- 
degenerate spectra, Attempts have been made to identify regions in the SUSY parameter 
space, which can accommodate the three patterns in turn. Our analysis clearly shows that 
the multi-dimensional parameter space of fiuSSM leaves enough room to accommodate all 
the diverse mass hierarchies of the three active light neutrinos. 

As a prerequisite of the loop calculation, we have derived the entire set of Feynman 
rules necessary for our analytical and numerical studies. The set of rules will be of immense 
importance, should one want to look for signatures of this model at the Large Hadron 
Collider (LHC), specifically through decays of the new particles of this model into SM 
particles. Approximate analytical forms of the entries of the expansion matrix have 
also been provided including all three generations of right handed neutrinos. A handful of 
relations within the four-component weak and mass eigenbasis for scalar and fermion fields 
of //z^SSM have also been worked out. 

The correlation of neutrino mass squared differences and mixing angles with relevant 
model parameters, consistent with the EWSB conditions, have been studied in detail for 
different mass hierarchies. The allowed regions of the parameter space are found to be 
rather seriously affected by the one-loop contributions. We have also observed that the 
globally fitted neutrino data can be accommodated into the /xi^SSM for various compo- 
sitions of the Lightest Supersymmetric Particle (LSP), once one starts playing with the 
parameters of the model. Although for the present study, we adhered to a Bino dominated 
LSP scheme, we numerically verified the possibility of having a light right-chiral neutrino 
(i.e. singlino) like LSP, compatible with neutrino data. This, by itself, is an interesting 
scenario in the sense that the singlino LSP offers a scope to probe the right-handed neu- 
trino mass-scale at the LHC. It would also be interesting to perform an explicit radiative 
correction to the heavy neutralinos. A dedicated analysis for this is beyond the scope of 
this paper, and a future publication with comprehensive discussion of all these issues may 
be well anticipated. In conclusion, one-loop radiative corrections to the neutrino masses 
and mixing angles for /xi^SSM are capable of substantially altering the tree level analysis. 
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A. Minimization equations 

The minimization equations with respect to the VEVs v'^, v^, V2, vi are given below. 

2J24'C' + E ^"'^c^i + J2i^ly'^j + + l^^'^2 + {A.xy = 0, (A.l) 

3 k j 

J2yJ^V2C^ + J^^mly'^^'j + E(^-^-)''^i^2 + Ig^vv'i + riv = 0, (A.2) 

3 3 3 

jit i 

-^X^V2C^ - H^riv'j - ^{Ax\yv1v2 + mj^^vi + -fg^^vi + fi'^vi = 0, (A.4) 

j 3 i 

where 

{A,xy = Y,{A,Y,y\]v2 + Y.iA^K^'^v^vl - {Ay}.yv^V2, 

3 j,k 

3 3 k 

i i 

V = Y r'A - fivi, p' = r' - X'vi, 

i 

l9 = \i9l+92), Cv = Yl^? + vl-vl 

(A.5) 

In deriving the above equations, it has been assumed that k^^^, {A/^ny^^, Yy-' , {AyYi,y^ , 
{m~cy^ , {rn-y^ are all symmetric in i,j,k. 
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B. Details of expansion matrix ^ 

In this appendix the entries of the expansion matrix ^ are given in details 



V2ginml.M2A ^ 

V2g2fimlcMiA 
^ UD 



6 



i3 



mlcM' 
2D 



I ^Xv2v'^ — 4:fj,A — + mycV2V%i — 3A [Xviv^ — 2mycv^V2)c 



. TJlycJ^f f. 2 ri o\2 21 

^ |Auiw Oj + m^^cwit; + 3A V2V Cij 



^^^^^ Uxf XvHl - -sm^2B) + '^v%hm2B + Avhm2B - AllMa] a 



2A I Au (1 sin 2/3) H 5-uS^sin2/3 + Av^sm2l3 - AfiMA ai - nmucv^cos2/3bi } , 



65 ~ 67 ~ ^jy- <X[Xv - AiJLMA)a2 H 62 - 2XfimycV2C2 

65 - 66 - {a (At;^ - 4/.M^)a3 + ^^^63 - 2A/.m..^;2c3} , (B.l) 



where 

ai = Y:'v2, hi = {Y^'vi + 2,Xv[), a = v[, 
•niyc = 2kv^, fi = 3Xv^, A = {kv'^'^ + XV1V2), 
V2 = vsin/3, vi = vcos(3, D = Det [M7X7] , 

M-M1 + M2' ^ "Hw-' ^^-^^ 

with i = e,/i,r =1,2,3. 

C. Scalar mass squared matrices 

The superpotential of the /ii^SSM violates Rp through lepton(L) number violation. This 
allows the Higgses (having zero lepton number) to mix with the sleptons (having non- 
zero lepton number). Hence, the neutral (both CP-odd and CP-even) and charged Higgs 
mass squared matrices are enlarged to 8 x 8, considering all three slepton generation. The 
independent entries of the CP-odd, CP-even and the charged scalar mass squared matrices 
were derived using eqs.(|A.3D, ([A.4p, and eq.(A.5). Details of each of these matrix elements 



were given in ref. [73], hence we do not repeat them here. However, we give the expressions 
for scalar quark (squark) mass squared matrices. 

In this appendix we present the relevant details of various scalar mass squared matrices 
required for the Feynman rules. The scalar sector of this model have also been addressed 
in ref. [72,73,79,80]. 
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C.l CP-odd neutral mass squared matrix 

In the weak interaction basis $p = {H^x^ -^2X' ^ni^ ^ni), the pseudoscalar mass term in the 
Lagrangian is of the form 

^psludoscalar = ^?M|,$p, (C.l) 

where Mp is an 8 x 8 symmetric matrix. The mass eigenvectors are defined as 

po = <;$p,, (C.2) 

with the diagonal mass matrix 

C.2 CP-even neutral mass squared matrix 

In the flavour basis or weak interaction basis = {H^-j^, H2-ji,u^-ji,i^n7i)-, the scalar mass 
term in the Lagrangian is of the form 

^TcXr = ^sMl^S, (C.4) 

where M| is an 8 x 8 symmetric matrix. The mass eigenvectors are 

5a = ^f>5,, (C.5) 

with the diagonal mass matrix 

iMt')lp = R'J,Ml,Rfs. (C.6) 
C.3 Charged scalar mass squared matrix 

In the weak basis <I>J = (if^^, i?^, ej^, ej^), basis the charged scalar mass term in the 
Lagrangian is of the form 

^charged scalar = ^C±^C^ (C-7) 

where is an 8 x 8 symmetric matrix. The mass eigenvectors are 

St = Rtp^%, (C.8) 

with the diagonal mass matrix 

One of the eight eigenvalues of the CP-odd scalar and the charged scalar mass squared 
matrix is zero and corresponds to the neutral and the charged Goldstone boson, respec- 
tively. We reiterate that explicit expressions for the independent entries of Mp, M| and 
are given in ref. [73]. 
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C.4 Scalar quark mass squared matrix 

In the weak basis, u'i = and d'i = {di^ , d*j^,) , we get 

^T,lZk = 1<^M'~ ^ + id^^Mj^ d'^ , (C.IO) 
where q = {u' , d'). Exphcitly for up and down type squarks (n, d), usmg eq.( [A.5 ) the entries 



are 



Q' 6' 2 2 



2 



n 

I 

(Mi)^^^^ = (Mi)^^^^ , (C.ll) 



and 



n 

n 

{Mj)RiLi = {M^)^i^' . (C.12) 
For the mass eigenstate q^ we have 

= Rlqj , (C.13) 

with the diagonal mass matrix 

{Mf% = RlMlR% . (C.14) 

D. Charged fermion mass matrix 
D.l Chargino mass matrix 

The lepton number (L) violation allows mixing between the MSSM charginos with the 
charged leptons and thus the chargino mass matrices enhances to 5 x 5. The mixing 
between MSSM charginos and charged leptons are governed by the left handed sneutrino 
VEVs(?;^) and neutrino Yukawa couplings. Both of which has to be small in order to satisfy 
the global neutrino data, thus this Rp violating mixing are very small. The chargino mass 
matrix for the /[/z/SSM have been addressed in ref. [72,73,79]. In the weak interaction basis 
defined by 

M/+^ = (-a+,^2^,e+,;U+,r+), 

*"^ = (-^A2 ,i?r,e^,;u^,r^), (D.l) 



n 2 
V 



1 ; 
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the charged fermion mass term in the Lagrangian is of the form 

/ 05x5 / 



r-mass 
^charged 



ym5x5 05x5 



(D.2) 



Here for simphcity we assume diagonal form of the charged Yukawa couphngs. The matrix 
"i5x5 is given by (using eq.( A.5[ )) 



/ M2 52^2 



"1-5x5 



\ 



g2Vi /i 











(D.3) 



V52< -rl Y^vi I 

The charged fermion masses are obtained by applying a bi-unitary transformation such 
that 

U^'m^^^V-^ =M%, (D.4) 

where C/* and V are two unitary matrices and M% is the diagonal matrix with non-negative 
entries corresponding to the physical fermion masses. The two-component mass eigenstates 
are defined by 



X~ = U,j^J, i,j = l,...,5. (D.5) 
Nevertheless, we notice that the 13, 14, and 15 elements of the chargino mass matrix (eq. 



(|D.3| )) are vanishing and given the orders of magnitude of various parameters, we also see 
that the values of the other off-diagonal entries (except for 12 and 21 elements) are very 
small. This indicates that the physical charged lepton eigenstates will have a very small 
admixture of charged higgsino and charged gaugino states. So we can very well assume 
(also verified numerically) that this mixing has very little effect on the mass eigenstates 
of the charged leptons. Thus, while writing down the neutrino mixing matrix, it will be 
justified to assume that one is working in the basis where the charged lepton mass matrix 
is already in the diagonal form. 

D.2 Quark mass matrix 

The mixing matrices for up and down quarks are 3x3 and they are diagonalized using 
bi-unitary transformation. Entries of up and down quark mass matrices ^13^3 and ^^^3 
are given below 



[m 



3x3Jij 



[m 



3x3)ij — "^l- 



(D.6) 
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The quark mass matrices are diagonalized as follows 



m 



3x3 



m 



u -1 
R 

d-1 



3x3-"-/? 



M 

M 



diag 
U ' 
diag 
D ■ 



(D.7) 



E. Feynman rules 

In this appendix we will study the relevant Feynman rules required for the calculations of 
the one- loop contributions to the neutralino masses. Some of the Feynman rules for this 
model have been derived in ref. [73]. Feynman rules for MSSM are given in ref. [113,114] and 
in ref. [115,116] for MSSM with singlet superfields. Feynman rules for iip-violating MSSM 
were studied in ref. [42]. The required Feynman rules are (using relations of appendix 
shown later) of the form neutralino-fermion- scalar /gauge boson and they are listed below. 

Neutralino- neutralino- neutral scalar 

The Lagrangian using four component spinor notation can be written as 



^mOltkPL + O^kPRW^Sl (E.l) 



where 



and 



(E.2) 



= {Olf^T- (E.3) 



Neutralino- neutralino- neutral pseudoscalar 

The Lagrangian using four component spinor notation can be written as 



^nna ^ ^q-^o^Pl + O^S^.Pr)^^ , (E.4) 



where 

nna 

Lijk 



Rf;(-^N:2N*3 + ^N*iN*3 



51 



V2 



V2 



and 



^Rijk — ~\yLijk) ■ 

Neutralino-neutralino-Z° 

The Lagrangian using four component spinor notation can be written as 



where 



1 



2 cose^ 



(Ni3N*3 - Ni4N*4 + Ni,^+7N* 

,m+7) ' 



1 



2 cose^^ 

Neutralino-chargino-charged scalar 

The Lagrangian using four component spinor notation can be written as 
where 



Lijk 



92 

+^^^2 ('^™U*2N^^^^4 — ^J""'U*^+2N|,n+4) 

+ •^fe.m+S ^52U*iN*^^_,_7 — — U*„j_,_2N*2 — --^\J*„^_^_2^ji 



9O 



ens 
Rijk 



^kl i^"^^i2^j,m+4 — ^ i,n+2i-'^j,m+7 

+kS (:^V,2N,-2 + -^V,2N,i +52VaN,-4 

+V^5lR-f,m+2Vj,m+2Nji 



2N, 



(E.5) 
(E.6) 



(E.7) 



(E.8) 



(E.9) 



(E.IO) 
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and 

Ones t r)cns \* r\ncs f/^cns \* /tti -i -i \ 

Lijk = yyRjik) ' ^Rijk = K^Ljik) ■ {^■^^) 

Neutralino-chargino-VF 

The Lagrangian using four component spinor notation can be written as 

= -g2xa^iO'EirPL + o'^^Pr)x',w+ - 92^,1^ [o^EiJPL + o-^JPr)xjW^. (E.12) 

where 



and 



OSJ = -U*iN,-2 - ^U|2N,-3 - -^U*„+2N,-,n+7, (E.13) 



Lij = V->Lji ) > ^Rij = V->Rji ) ■ (E.14j 



The factors rjj and are the proper signs of neutrahno and chargino masses [115]. They 
have values ±1. 



Neutralino-quark-squark 

The Lagrangian using four component spinor notation can be written as 

C-'^ = -m{Ol%PL + 0^^lPn)x',qk - mO'nJkPL + 0^^lPn)qjrk. (E.15) 

where 

0'ESk = iO'g^*^ 0^gk = (0^kr, (E.16) 



and 



~/^nuu \/-mn*-riu -ivt tju* tvt pu* 

9'^Rijk — ^km^^i4:-^Rjn ~ '^^^k,m+3^^il-^Rjm^ 

—r^ndd -r>d I Jl^-yj* "od , TM"* D<i I _L vnm-r)d tvt* "od 

9'-^Lijk — "-km \~~^^^i2"'Ljrn ~^ ^^^^il"^Ljm J ~^ ^d "^k,m+3^^i3"'Ljn> 

~/^ndd vmn*T}d xT 7?<^* i ^ffl T}d tvt E?<^* ^"1? 1 

a^Rijk- ^km'^iS^Rjn + '^;y^^k,m+3^^il^Rjm- 1^.17) 
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F. The and n]< function 

In this appendix we give the detail expressions for the renormahzed self energy functions 
TY, and n,^. The net result is 



167r2 



~2 8 10 



r=l k=l 



m~o,mso) 



7^ 



7 10 



y \ ^ \ ^ f /^nna r\nna i r\nna rynna \ jd ( J2 ™2 _,2 \ 
~ Y 2^ 2^ y^Lkir^Rjkr + ^ Ljkr'-' Rkir) ^nP > ™^ > "T-poj 
r=l k=l 
10 

+ 92 2^ \'->Lki ^Ljk + ^Rki ^Rjk) ^1 KP > ^yll > ) 
k=l 

k=l ^ 



+ 



r-=l fe=l 

E E [oT^irOT^kr + Ol%0^l) B,{p',ml,ml) 

r=l k=l 
6 3 



r=l k=l 



+ ^9^ E E ( ^^fci*- j fcr + O^Ljkr Ofofir ) Bi{p^,m'^^, 171- 



(F.l) 



n, 



167r2 



8 10 



r=l A;=l 
7 10 



nir-o 



~ 9 Z^Z^^^[^Lkir^Ljkr + ^Rkir^Rjkr)Bo{P ,m~o,mpo) 
r=l k=l 
10 

o _2 \ ^ {^nnz^nnz , /-\nnz^nnz\ tj ( Jl ^2 ^2 \ 

- ^52 2^ "^xg l<^Lfei '->Rjk + <^Lifc <-'Rfei j ^0 , rn~o^ , m^o ) 



k=l 

5 



fc=i 
7 5 



r=l fc=l 

6 3 

^3 2^1^ rriuk [yLkir^Ljkr + ^Rkir^Rjkr ) ^OKP > 



r=l k=l 
6 3 



r=l k=l 



(F.2) 



Detail expressions for the couplings are given in appendix The Bq, Bi functions are 
given in appendix The factor 3 appearing in front of the quark-squark loop contributions 
signifies three variations of quark colour. 

G. The Bq and Bi function 

The Bq and Bi functions are Passarino-Veltman [94, 95] functions defined in the notation 
of [96] 

^ r> ( 2 2 2 \ _ 4-D f 1 \ 

^^^^l^^[p ,mf,,m,j - J (2vr)0(^2_^2 )((^+^)2_^2)' 

Bf, {p^ , ml, ,ml ) = p^Bi {p^ , ml, ,ml). 

(G.l) 

H. Some Useful Relations 
Fermionic sector 

For neutralinos the following relations between mass and weak eigenstates are very useful 



PL^k = PLN*k+7xl Pli^I = PLN*k+4xl 

Pr& = PnNiiXi, PrW^ = PRNi2Xi, PrHj = PRNij+2, 

PRVk = PRNi^k+ixt Pr^I = PRN^,k+ixl 

where j = l,2 and A; = 1,2,3, (H.l) 



with i varies from 1 to 10 and 



^L=(^V PR=m^]- (H.2) 



In terms of the four component spinors Xi for charginos, the following relations between 
mass and weak eigenstates are very useful. 



PlW = PlV*^x^, 

PrW = PRUilX^, 

PlW'' = PlU:,xI, 
PrW' = PrVuxI, 



PlH = PLV*2X^, 
PrH = PRU,2X^, 

PlH"" = PlU^xI, 

PrH' = PRV,2Xt, 



^dk = PLUlk+2Xi, 
PRh = PRV^,k+2Xt, 
Pdk = PLV*k+2Xi, 
Pr^I = PRUi,k+2Xii 
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where k = 1,2,3, and i varies from 1 to 5. The last six relations are for the charge- 
conjugated fields. 

The four component neutralino, chargino and charge conjugated chargino spinors are 
respectively defined as 



(H.4) 



where and xf are two component neutral and charged spinors, respectively. 
Scalar sector 

The relations between weak and mass eigenstates for neutral scalar, neutral pseudoscalar 
and charged scalar are given by eqs.( p.2| ), ( |C.5 ), and ( |C.^ ). 
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